Implication of possible observation of enhanced /jl^/jl decay 
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The very rare — >> decay may be the last chance for New Physics in flavor sector at 

the LHC, before the 13 TeV run in 2015. Partially motivated by the known tension in sin2/3/0i, 
enhancement beyond (3-4)xl0 -10 would likely imply the effect of a fourth generation of quarks. 
If observed at this level, the 126 GeV boson may not be the actual Higgs boson, while the b — »• d 
quadrangle (modulo m t /) would jump out. The 2011-2012 data is likely not sensitive to values below 
3xl0 -10 , and the mode should continue to be pursued with the 13 TeV run. 
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I. INTRODUCTION 

Despite the discovery |l|, [2| of a 126 GeV boson by the 
Large Hadron Collider (LHC) in 2012, the LHC has so 
far been a disappointment: no New Physics beyond the 
Standard Model (SM) has been seen, and even the new 
boson appears Higgs-like, i.e. as prescribed by SM. 

Surveying the terrain, there seems one last hope for 
discovering New Physics, namely — ^ There 
is some motivation for enhancement, from the well 
known [3|, |4[ mild (of order 2<j) but lingering tension be- 
tween direct measurement of CP violation (CPV) phase 
of Bd-Bd mixing, versus extraction by indirect means. If 
an enhanced — »> /i + \i~ rate is discovered with 2011- 
2012 LHC data, the likely explanation would be a fourth 
generation of quarks. This would then cast doubt on the 
Higgs boson interpretation of the 126 GeV boson. 

The B® — >> decay has been a highlight pursuit 

since Tevatron times, and only recently surpassed [5| in 
sensitivity by the LHC. The drive has been the possibly 
huge enhancement by exotic scalar effects inspired by su- 
persymmetry (SUSY), but now excluded by the first ev- 
idence for SM-like rates by the LHCb experiment [6| . In 
contrast, the search for B^ — >> /i + \i~ has not shared the 
limelight. This is because the SM prediction itself is 30 
times lower than B® — >> However, the combined 

LHC bound is now within [B] a factor of 8 of the SM 
prediction, and one may ask whether this mode could be 
anywhere enhanced up to this order. 

As pictorialized by the "Straub plot" [7[ and discussed 
recently by Stone (8), most models of enhancement for 
B® — ^ have now been eliminated by the SM- 

like B® — ^ rate measured by LHCb, with two 

exceptions. One is an old, purely left-handed SUSY 
model [9]. However, the region allowed by current data 
is but a corner of the parameter space, hence not plau- 
sible. The other would be [10] the 4th generation (4G), 
where B^ — >• /i + \i~ and B® — >> decays are mod- 

ulated by different Cabibbo-Kobayashi-Maskawa (CKM) 
products V t 1 d Vt'b and V£ s Vt>b, allowing B® — > ji + ji~ to be 
enhanced up to the current bound, even if B® — ^ is 
SM-like. Stone has followed conventional wisdom to ar- 
gue [8[ that 4G has been "eliminated by the Higgs discov- 
ery" , because it "would cause the Higgs production cross- 



section to be nine times larger ..." [ll|. In fact, a com- 
prehensive analysis [12| including electroweak and flavor 
observables plus earlier Higgs production data already 
ruled out 4G in SM framework. There are two catches 
in this pessimism, however. First of all, it is not yet es- 
tablished that the observed 126 GeV object is the Higgs 
boson of SM. For example, a dilaton might mimic [13j 
the Higgs with current data. Second, the Higgs boson 
of SM does not enter into the B® — ^ /i+ fi~ process (the 
same holds for the B^ box diagram and B + — ^ 7r + /i + /i~ 
processes we consider). To assume indirect arguments in 
the flavor pursuit is self-defeating, especially when there 
is still room for large enhancement; it actually highlights 
the potential impact of a discovery. 

It was shown [14[ recently, through an empirical gap 
equation [l5], that dynamical electroweak symmetry 
breaking (DEWSB) could occur through strong Yukawa 
coupling of 4G quarks. Although there is no account 
for how a dilaton actually emerges, the scale invar i- 
ance of this gap equation allows for a dilaton to appear. 
The dilaton possibility can be checked experimentally 
through the absence, or suppression, of vector boson fu- 
sion (VBF) and associated production (VH) processes, 
which requires more data than currently available. The 
very large Yukawa coupling needed for DEWSB is consis- 
tent with not finding the 4G quarks so far, where the cur- 
rent bounds [l6| are already above the nominal [ItJ uni- 
tarity bound (UB). Thus, the numerical study we present 
below is only meant as an illustration. 

In the following, we review input parameters and con- 
straints, then present our numerical study. We indeed 
find enhancement beyond 4 x 10 -10 (4 times SM) is pos- 
sible [18] within the parameter space indicated by the 
known tension in sin2<I>^ d = sin20i//3. We give an as- 
sessment of immediate and longer term prospects. 



II. CONSTRAINTS AND INPUT PARAMETERS 

There is no indication for New Physics in b — »• s tran- 
sitions at present. The best probe is sin2$^ s measure- 
ment pursued by LHCb, where <£# s is defined as the 
CPV phase in the B s — >> B s mixing amplitude (hence 
sin2<l>B s = sin S ). This definition is consistent with 
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sin 2/3/01 = sin 2$^ used by the B factories. The 4G 
t' quark could have easily affected many b — >• s pro- 
cesses [13, However, all of these, including s — d 
transition effects, can be tuned away or softened by a 
small iV^^/frl strength, which is demanded by sin2<3>£ s 
being consistent with SM expectations and is yet to be 
measured. As illustrated by the Straub plot [7|, llOj , 
B s — )> and Bd —> can vary independently 

from each other, i.e. through V t * d Vt'b and V t * s Vfb: sub- 
ject to constraint from kaon physics (affected by V t 1 d V t ' s ). 

It is well known [3|, |4[, however, that there is some 
tension between the directly measured value [20[ of 

sin 2/3/01 = 0.679 ±0.020, (1) 

and SM expectation via p/fa = argA^ M , where [H| 



A 



SM 



-K -K- -\V ud \\V ub \e-^ + \V cd \\V cb l (2) 



with Xi = V* d Vib. The terms on right-hand side of Eq. (2) 
can be measured at the tree level. Currently [20j , 



(68 



10) 



(3) 



and we take the central values \V u d\ = 0.974, \V cd \ =0.23 
and \V c b\ = 0.041 [20|. Variations in these values are not 
central to our discussion. 

In contrast, \V u b\ also has some tension in the measured 
values. Extraction via inclusive or exclusive semilep- 
tonic B decavs yield approximately 4.41 x 10 -3 and 
3.23 x 10 -3 [20], respectively, with the average value of 
4.15 x 10 -3 (the inclusive approach has better statistics). 
We use central values, as our purpose is only for illustra- 
tion, hence we will treat the average (which is close to 
inclusive) and exclusive cases separately. 

Although the strength of |A £ SM | ~ 0.0088 is not sensi- 
tive to \V u b\, the phase is sensitive to its value, 



sin 2/3/</>i = 



0.76 for |Kb| ave 
0.63 for|Kfcl excl , 



(4) 



which both deviate from Eq. (1) by more than 2a (the 
inclusive value of 0.81 deviates even more). This devi- 
ation offers some motivation for New Physics in b — >• d 
transitions. It could easily be due to the 4G quark t' ', 
where one simply augments Eq. (2) by 



At — A^ M — \f , 



(5) 



and the b — > d triangle becomes a quadrangle 



A n + A c + A, + A t / = 0. (6) 
In our following study, we parameterize [22| 



A*' = r db e % 



(7) 



In our phase convention, A c = V* d V cb is practically real, 
while A n = V* d Vub is basically the same as in SM. 

To study sin2<I>£ d and B(B d — >• /i + /i~) in the Tdb-^db 
plane, other constraints should be considered: 



• radiative b — d^y processes (including B — ^ pj) 
is ineffective because it is hard to separate from 
b — >• 57, difficult to study with LHCb, and in any 
case insensitive to virtual 4G effects; 

• B — )> 7T7T decays, while quite well studied, suffers 
from hadronic effects (even B — >> Ktt suffers from 
hadronic effects), and do not provide good con- 
straints; 

• the well measured ArriB d provides a constraint 
through uncertainties in f% BB d \ 

• only very recently was the electroweak penguin 
B + — )> 7r + /i + /i~ decay measured [23[, in contrast 
to electroweak b — > s penguins. 

Although it may be a little surprising, there are not many 
observables that provide sound constraints on \ t r. We 
collect below the relevant formulas for our study. 
The t' effect to Bd mixing 



Am Bd ^ 
sin2$£ d 



G%M ™ m Bd B Bd f Bd VB \Af 2 \, 



6tt 2 

: sin(argAf 2 ), 
is (explicit forms can be found in Ref. pi} ) 



(8) 



A d — 
^12 — 



(A t SM ) 2 5 (x t ) 



2X^ M X t >AS^ 



AS, 
AS, 



(i) 


(2) 



where Xi = m 2 /Mw 



\ 2 t ,AS£\ (9) 
So(xt,xt>)-So(x t ), (10) 
So(xt>)-2S (xt,xt>) + So(xt), (11) 
Besides 4G parameters, the main 



uncertainty is in [25] 



/ Bd B^ 2 = (227±19) MeV. 



(12) 



For the current bound [5[ of 

B(B d < 8-1 x 10" 



-10 



(13) 

our purpose is to illustrate whether, and how, it could 
get enhanced to such values by 4G effect. Here, we use 
the usual trick |26j of "normalizing" the branching ratio, 



B(B d ->• M +/0 



_ B{B d -»• n+ir) 



c 



Am Bd 

T Bd Am e *v ft \\f M Y (x t ) + \ t ,AY \ 2 



(14) 



where A>o = 
Ref. [10], and 



B Bd Vb |Af 2 | 
= Y (x t >) - Y (x t ) with Y (x) given in 



C = 6tt 



47r shr 6 



w 



(15) 



Through the ratio of Eq. (14), one not only eliminates the 
hadronic parameter fs d , but the A^ M factor also cancels 
in the SM case, and one recovers the SM result of 1.1 x 
10 -10 , with little sensitivity to \V ub \. 

The treatment of B + — >• 7r + /i + /i~ would be given in 
the next section. 
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<f>db = arg(V t *, d Vf b ) 4> db = arg(V t *, d Vf b ) 

FIG. 1. Allowed region in |V^ d Vi/b|-arg V t */ d V t 'b (i-e. rdb-<frdb) plane for (a) average (b) exclusive \V u b\ values, for m t ' — 700 
GeV. The solid-blue lines are labeled 10 10 B(Bd — > contours, where above the value of 8 (semi-transparent gray) is 

excluded by the combined result of LHC experiments. The dark (light) narrow green-shaded contours correspond to the l(2)a 
regions of sin (Eq. (1)), while the broad pink-shaded contours correspond to the l(2)cr regions of ArriB d allowed by 
Eq. (12). 




III. PHENOMENOLOGICAL STUDY WITH 
HEAVY t' 

We plot in Fig. 1 for my = 700 GeV the 2a range in the 
Tdb-^db plane, for sin 2&B d (green) allowed [27] by exper- 
imental measurement of Eq. (1), Am# d (pink) allowed by 
lattice error in Eq. (12), and the bound on Bd —> 
(gray exclusion) according to Eq. (13). We include la- 
beled contours of 0.5, 1, 2, 4, 8 for 10 10 B(B d /i + /i"). 
Fig. 1(a) and (b) are for taking \V u b\ to be the central 
values of 4.15 x 10" 3 and 3.23 x 10~ 3 , respectively, for 
the mean (between inclusive and exclusive) and exclusive 
values from semileptonic B decay studies. 

Consider Fig. 1(a), i.e. for \V u b\ = 4.15 x 10 -3 , the av- 
erage between inclusive and exclusive measurements (the 
inclusive case is qualitatively similar). The well measured 
CP phase sin is sensitive to t' effects, but free from 
hadronic uncertainties, hence the narrow (green) contour 
bands. In contrast, Am# d is less sensitive to <pdb, and 
more accommodating because of hadronic uncertainty in 
fs d B B ^ . The broad (pink) contour bands show the 1 and 
2<T allowed region by Eq. (12), and rules out a branch of 
the sin2<3># d contour (for <pdb between —10° to 15°), due 



to coherent enhancement of Am^ from t' effects. 

Consider now the gray excluded region from the com- 
bined LHC bound on Bd — >• Eq. (13). It is seen 
that there are two slivers of parameter space, around 
(Ub, <t>db) ~ (0.0025, 180°) (region A) and (0.002, 252°) 
(region B), where B(Bd — » M + aO could be above 4 x 
10 -10 , or enhanced by 4 times over SM, which are dis- 
covery zones for 2011-2012 LHC data. Near region B, 
B(Bd — » quickly drops below 4 x 10 -10 as Tdb 
becomes weaker than 0.002. For (j>db ~ 245° and Tdb 
varying from 0.0008 to 0.0015, B(B d -> hovers 
at (1-2) xlO -10 , while for r db ~ 0.0004 to 0.0008 and 
(j)db varying from 240° to 330°, B(Bd —> hovers 
at (0.5-2) xlO" 10 , i.e. within a factor of two of SM ex- 
pectations. These regions, combining to a broad crescent 
shape which we refer to as "region C" , would likely need 
much more data to probe. 

The LHCb experiment has recently measured [23j] 

B(B + 7T = (2.3 ± 0.6 ± 0.1) x 10~ 8 , (16) 

which is the rarest B decay observed to date. The result 
is consistent with SM expectations, but interpretation 
depends on form factor models. To reduce form factor 
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FIG. 3. Same as Fig. 1, but for m t > = 1000 GeV. 




4>db = arg(V t *, d V t < b ) 4> db = arg(F; rf F f . fi ) 

FIG. 4. Same as Fig. 2, but for m t , = 1000 GeV. 



dependence, we take the ratio 



Bttl 



B(B+ 



4G 



£(£+^7T+/i+/i-)| 



SM 



(17) 



where both 4G and SM results are integrated from 
q 2 = (1. 6) GeV 2 , which is under better numerical con- 
trol [23, [29j. Since this does not match what LHCb 
does, we draw contours in Fig. 2 (red-dashed), and view 
R-Kfifi ~ 2-3 as the range beyond which LHCb would have 
found inconsistency with SM expectations. Thus, we are 
interpreting LHCb's statement of consistency with SM, 
allowing for form factor uncertainties. It is clear that 
this approach is not as good as the zero crossing point 
for Afb(B —> K*jjLjjL), but this is the first observation of 
rare b — > dlt decays, compared to the decade- long explo- 
ration of b — >• sit processes. For numerics, we combine 
Wilson coefficients at next-to-leading order with leading 
order decay amplitude based on the QCD factorization 
approach [28|, [29j. For dealing with New Physics, and as 
we take a ratio, this should suffice for our purpose. 

If we now compared Fig. 1(a) with Fig. 2(a), we see 
that ArriB d is more powerful than B(B + — ^ 7r + /i + /i~) 
in excluding the sin 2 §B d -allowed branch near <pdb ~ 0. 
This is reasonable, since B + — >> 7r + /i + /i~ is only recently 
observed and prone to hadronic form factor uncertain- 
ties, while Ara# d has been measured since 25 years, 

~ 1 /2 

with hadronic uncertainty narrowed down to fs d B^ , 
which itself has been subject to intense lattice studies for 



years. It is, however, comforting to see that for region 
A, R-Kfifi is not more than 2 (except the upper reach near 
<pdb ~ 190°), hence should be easy to accommodate by 
form factors, while for regions B and especially region C, 
R-Kfifi is even less than 2 and closer to 1. Thus, the newly 
measured B + — » 7r + /i + /i _ does provide a sanity check. 

Turning to the case of exclusive \V u b\ value, Fig. 1(b) 
and 2(b), we find that regions A and B basically switch 
roles. This is because for \V u b\ ~ 3.23 x 10 -3 , the 
expected sin 2$^ value in SM falls below that of di- 
rect measurement, as seen in comparing Eq. (4) to 
Eq. (1). Calling it region A', the sliver of region around 
(rdb, <t>db) ~ (0.002, 160°) could enhance B(B d H> /i + /i") 
more than 4 times above SM, and observable with present 
LHC data. Region A' extends to the broad crescent re- 
gion C, where even rdb values as lower as 0.0002 could 
account for the measured sin2$B d , but B(Bd —> /i + /i~) 
can be probed only beyond 2015. Again, ArriB d ex- 
cludes the sin 2$ s d - allowed branch around <fidb ~ 30°. 
Region B' is now a considerably broader region in param- 
eter space that allows enhancement of B(Bd — » 
above 4 x 10 -10 . For example, for rdb above 0.0023 and 
<pdb above 230°, B(Bd — » /i + /i _ ) can be greater than 
6 x lO -10 , fB d B B l2 d is within 2a of Eq. (12), while 
is not more than 2. We also see that, for region B', R nfIfI 
provides as good, perhaps better constraint, than Ara# d , 
disfavoring the region of rdb greater than 0.0025 around 
4>db ~ 205°, that seems perfectly allowed by Am# d . 

Now let us consider mf values. The 700 GeV value 
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FIG. 5. Sample b -> d quadrangles for A t / = V t 1 d V t , b = 0.0025 e il80 ° with average \V ub \ = 4.15 x 10" 3 (left), and for 
\ t , = V*, d V t , b = 0.0023 e i230 ° with exclusive \V ub \ = 3.23 x 10" 3 (right). 



used so far is just above current experimental limits [l6j . 
and correspond to Yukawa coupling strength y t > ~ 4, or 
af — 1.3, which is why there is UB violation (UBV). 
However, we do not quite know what is the true ex- 
pansion parameter. Furthermore, even if perturbation 
breaks down, it does not mean there is no t f effect. In 
fact, perturbation in \ t > certainly holds, though the func- 
tions ASq^ and AYb in Eqs. (9) and (14) gets modified by 
UBV effects. The overall form of these equations should 
not change. We therefore consider the r?v — 1000 GeV 
case, i.e. a t > ~ 2.6, to illustrate the situation far beyond 
UBV [17]. Note that Ref. [14] finds DEWSB occurs for 
Dq (the 4G doublet is treated as very close to degenerate) 
of order 47r, i.e. of order the ttNN coupling, implying 4G 
quark masses no less than 2 TeV! 

The plots corresponding to Figs. 1 and 2, but with 
7?V = 1000 GeV, are given in Figs. 3 and 4. We gener- 
ally see reduced rdb values. Region A is now excluded, 
but regions B, A' and B' become more robust in Am^, 
but values for B(B d -> /i + /i~) higher than (5-6) xlO -10 
are slightly disfavored by B + — » 7r + /i + /i~. Viewed differ- 
ently, if enhanced Bd — > is discovered, one may try 
to scrutinize whether B + — > 7r + /i + /i _ is also somewhat 
enhanced beyond SM. Regions C and C generally stand 
well, with at best mildly enhanced Bd — ^ 



IV. DISCUSSION AND CONCLUSION 

We started with the question of what could still en- 
hance Bd —> decay, when everything at the LHC 
seems consistent with SM. The answer is that, proba- 
bly only the 4G t f quark could do the job, even if 4G 
seems disfavored by the Higgs-like nature of the 126 GeV 
boson. Admittedly, even if 4G is the explanation for 
the sin2$^ d tension as seen by the B factories, to have 
B(Bd — » ji + to be within a factor of 2 of the current 
bound of 8.1 x 10 -10 is only a fraction of the allowed 
parameter space, hence not particularly likely. However, 
only with such enhancement is there any chance for LHC 
experiments to make the discovery with 2011-2012 data, 
and discovery it indeed will be. If discovered — within 
2013 — then not only 4G would get uplifted, some doubt 
would be cast on the SM Higgs nature of the 126 GeV 
boson, while "impostors" such as dilaton would gain in 
weight. We have remarked in the Introduction that it 



would take the establishment of VBF and VH produc- 
tion processes to exclude the dilaton possibility, which 
cannot be achieved with 2011-2012 data [14j . 

An intriguing outcome of discovering Bd — » de- 
cay would be that, all of a sudden, the b — >> d triangle 
falls into our lap! Let us illustrate. Since r?v = 1000 
GeV cases have smaller Tdb = \Xf\ = l^d^t'fcl values, for 
reasons of plotting, we take two examples from r?v = 700 
GeV. From region A of Fig. 1(a) (average \y u \\ = 4.15 x 
10- 3 ), we take X t > = V t 1 d V t > b = 0.0025 e il80 °. From re- 
gion B of Fig. 1(b) (exclusive \V ub \ = 3.23 x 10~ 3 ), we 
take X t > = V* d V t > h = 0.0023 e i230 °. 

The quadrangle of Eq. (6) is constructed as follows. 
To simplify discussions, we normalize to A c = V* d V c b = 
—0.0094, which becomes a unit vector pointing left. 
Then, A n = V* d V uh /\\ c \ = 0.44 e - i68 °, 0.34 e -* 68 °, re- 
spectively, for the average and inclusive cases, with cor- 
responding A t / = 0.27e il80 °, 0.24e i230 °. Then X t just 
connects the tip of A n with the end of Xf . The two ex- 
amples for 700 GeV are plotted in Fig. [5] in the form 
to compare with the usual SM triangle [20[. These are 
relatively precise quadrangles, and illustrate how 4G ac- 
counts for a shift in sin2$^ d away from SM expectation, 
where is the angle between the dashed line, A^ M and 
the real axis. Since t' is much heavier than t, a smaller 
\f could cause the shift. 

The sample b — >• d quadrangles are for largest allowed 
solutions for 7^5, i.e. regions A (for |ln6| ave ) and B' (for 
|U n6 | excl ) for mt> = 700 GeV, and would be the case if 
Bd —> is discovered soon. They are relatively ex- 

treme, however, since even for m t > = 700 GeV, regions 
C and C can provide solutions for sin for much 

smaller rdb = \V t 1 d Vt>b\ values, with possible phase values 
extending over a large range. For heavier t' illustrated by 
1000 GeV, \V^ d V t >b\ is smaller by half compared to 700 
GeV case, with region A is eliminated. 

The quadrangles of Fig. [5] reminds us of the possi- 
ble [30[ link to the baryon asymmetry of the Universe 
(BAU): 4G greatly enhances CPV from SM, and is seem- 
ingly sufficient for BAU (although a first order phase 
transition remains an issue), which boosts the merit of 
4G. It does not depend much on the area of the quadran- 
gle, as the enhancement rests in powers of r?v and . 
We note that X t f in Fig. [5j though smaller in strength 
than At and A c , is not that small compared with X u . 
Furthermore, we know that \Vfb\ cannot be more than 
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0.1 [3l|, especially for our large m t ' values. Hence, |Af| 
plotted in Fig. [5] correspond to \Vfd\ that is larger than 
\V t d\ SM ^ 0.0088, which does not fit the CKM pattern of 
trickling off as one goes further off-diagonal. One could 
use this to argue that enhanced Bd — > decay to 

the level observable with 2011-2012 data is not plausible. 
However, the issue is best left to experiment. 

For m t > = 1000 GeV, \\ t ' | values tend to drop by half, 
but \Vfd\ would still be comparable to \V t d\- Only if 
one gives up enhancement would the ratio \Vt'd/Vtd\ turn 
"natural". In fact, for the exclusive value case for V u b, 
| A*/ 1 (i.e. rdb) could be (1-2) xlO -4 and still account for 
sin2$B d "anomaly". Such values for \Vfd\ would become 
"natural" when compared with \Vtd\- However, even if 
4G gains support by 2015, this region (C and C) would 
need a very large data set to explore. 

We conclude that 2013 remains a pivotal year where 
one could discover the very rare Bd — » decay 
mode at over 4 times SM expectations. The chance is 



not large, but not zero either, with partial motivation 
from the (mild) sm2&B d discrepancy. If discovered with 
2011-2012 data set, the implications would be quite huge: 
uplifting the 4th generation (with prospect of CPV for 
BAU), casting some doubt on the SM Higgs interpreta- 
tion of the 126 GeV boson, and perhaps the only New 
Physics (at least in flavor sector) uncovered at the 7 and 8 
TeV runs at the LHC. But it is more likely that the LHC 
would once again push the limits down towards SM. If 
such is the case, the fate of the 4G would have to be de- 
termined elsewhere. But Bd — >• should certainly 
be pursued further at the 13 TeV run. 
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